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We present a multiscale quantum-defect theory based on the first analytic solution for a two-scale 
long range potential consisting of a Coulomb potential and a polarization potential. In its application 
to atomic structure, the theory extends the systematic understanding of atomic Rydberg states, as 
afforded by the standard single-scale quantum-defect theory, to a much greater range of energies 
to include the first few excited states and even the ground state. Such a level of understanding 
has important implications not only on atomic structure, but also on the electronic structure of 
molecules and on atomic and molecular interactions and reactions. We demonstrate the theory by 
showing that it provides an analytic description of the energy variations of the standard Coulomb 
quantum defects for alkali-metal atoms. 

PACS numbers: 31.10.+z,32.80.Ee,32.30.-r,31.90.+s 


The Rydberg-Ritz formula jj], [§] for atomic spectrum 
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was one of the very first universal properties uncovered 
for quantum systems, and played an important role in 
the very establishment of the quantum theory. Here E- mn 
represents the ionization energy, Rm = Roo/ (1 + m e /M) 
is the reduced mass Rydberg constant (with M being 
the ion mass), n is the principle quantum number, l 
and j are the orbital and total angular momenta of the 
electron, respectively, and /r^ oul is the Coulomb quan¬ 
tum defect. The formula asserts that despite consider¬ 
able differences in atomic spectra, the Rydberg series for 
different atoms differ from each other only in a quan¬ 
tum defect which encapsulates all complexities of short- 
range interactions. This universality, which originates 
from the fact that a sufficiently highly excited electron 
sees mostly the Coulomb potential, stimulated the devel¬ 
opment of the quantum-defect theory (QDT) and mul¬ 
tichannel quantum-defect theory (MQDT) 0-0]. They 
have long become the standard for understanding atomic 
and molecular spectra and electron-ion interactions 
The universality as represented by the Rydberg-Ritz 
formula with a constant il g ° u1 is however strictly appli¬ 
cable only to sufficiently highly excited Rydberg states. 
This is reflected, especially for atoms with highly polar¬ 
izable cores, by a significant energy dependence of /ry oul 
for lower lying states (see, e.g., Refs. [8j-t.ll|)- Does there 
exist a more general universality that applies also to the 
first few excited states or even the ground state, which 
are often of more practical and experimental interest? 
The answer to this question has implications far beyond 



atomic structure. Not only will it determine the degree 
to which other single-atom properties, such as oscillator 
strengths, ionization cross sections, static and dynamic 
polarizabilities (see, e.g., Ref. [ill), follow universal be¬ 
haviors outside of the Rydberg regime. It will also deter¬ 
mine the degree to which important atomic interaction 
parameters, such as the C$ van der Waals coefficients 
(see, e.g., Ref. 0 ) , follow universal behaviors for differ¬ 
ent atoms and different electronic states. Even further, 
it will determine the relations among different electronic 
states of a molecule and the relations among electronic 
states of different molecules. The prospect for such sys¬ 
tematic understanding of an entire manifold of molecular 
electronic states can be crucial for understanding atomic 
and molecular interactions and reactions, especially in 
excited electronic states where many of them participate 
simultaneously. 

This work establishes a broader universality in atomic 
structure as the first application of an analytic two- 
scale QDT. The theory is based on our analytic solution 
for the “Coulomb+Polarization” potential of the form 
—C\/r — Ci/r 4 . It is, to the best of our knowledge, 
the first analytic solution of the Schodinger equation for 
a two-scale central potential, for which no analytic so¬ 
lution is previously known or expected. Our solution 
represents a new class of special functions and has po¬ 
tential for generalization to other multiscale potentials. 
In applying the two-scale QDT to atomic structure, we 
show that the theory introduces a “new” quantum defect 
that has much weaker energy dependence than the tradi¬ 
tional Coulomb quantum defects fj, Goul that goes into the 
Rydberg-Ritz formula, and the theory provides an ana¬ 
lytic description of the energy dependence of /j, Coul for Rb 
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and Cs atoms, down to their ground states. The results 
establishes, in an analytic framework, a broader univer¬ 
sality in atomic structure and spectrum. Physically, it 
means that the statement that a Rydberg electron sees 
mostly the Coulomb potential can be replaced by a more 
general statement that an excited or an outer electron 
sees mostly a Coulomb potential plus a polarization po¬ 
tential, to a remarkable accuracy. 

We note that the importance of polarization in atomic 
structure is well known, as reflected both in perturba¬ 
tion calculations for alkali metals [13, EH j and in model 
potentials chosen, e.g., both for “one-electron” alkali 
atoms [3E3 and for “two-electron” alkaline earth atoms 
i E3- The underlying universality was however diffi¬ 
cult to identify, define, or describe, since a perturbative 
treatment is only applicable to high angular momentum 
states, and the effect of core polarization is difficult to 
distinguish from other short-range effects in a numerical 
calculation. 

Our two-scale QDT for “Coulomb+Polarization” po¬ 
tential is built upon the analytic solution of the radial 
Schrodinger equation 
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where /i is the reduced mass. It has two characteristic 
length scales p n := (2 ^,C n /h 2 ) 1 ^ n ~ 2 ' ) (n = 1,4), corre¬ 
sponding to each one of the potential terms of the form 
of —C n /r n (n = 1,4). Each length scale /?„ has a corre¬ 
sponding energy scale of 4 B) = ( fi2 / 2 M)(l//3n) 2 - Scaling 
the radius r by Pi, and the energy e by its corresponding 
Sg\ the scaled Schrodinger equation takes the form of 
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dr 2 
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where r s := r/P\ and e s := e/sg \ with the ratio of length 
scales /3 4 //3\ serving as a measure of the relative strength 
of the Coulomb and the polarization potentials. 

Through generalizations of techniques that have led to 
other single-scale QDT solutions for l/r 6 jl9j], 1 /r 3 [20j |. 
and l/r 4 (2ll - [25| potentials, we have solved Eq. © an¬ 
alytically to obtain its QDT base pair of solutions and 
the corresponding QDT functions [26[. The base pair /° 
and g c are defined with energy and partial wave indepen¬ 
dent asymptotic behaviors around the origin (specifically 
r s < Pi/Pi), 


fliirs) ~r s J cos [(Pi/Pi)/r s - tt/4] , (4) 


9ej( r s) - r s ^j sin KM Pi)/ r a - *V 4 ] . (5) 

for all energies e s [2|| [23] . They are normalized 
such that their Wronskian W(f c ,g c ) := f c (dg c /dr s ) — 


( df c /dr s )g c = 2/7r. For e s < 0, we have shown that the 
QDT base pair has asymptotic behavior at large r s given 

by 
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where k s = (—Cs) 1 / 2 . It gives a 2 x 2 W c matrix with 
elements W£ y describing the evolution of a wave function 
through the — C\/r—C 4 /r 4 potential at negative energies 

M- 

In terms of the W c elements, the bound spectrum 
for any potential V (r) that behaves asymptotically as 
—Ci/r — C 4 /r A can be formulated f2(| as the solutions of 


Xi( e s,Pi/Pi) = K c (e,l,j) . (8) 


Here = W(/_/Wg_ is a universal function of the scaled 
energy e s and the ratio of length scales P 4 /Pi- The K c is 
a short-range K matrix defined by matching the short- 
range wave function, u e i(r) for potential V(r), to a linear 
combination of the QDT base pair [26[ [27j 


Ueij(r) = A elj lfh(r s ) - K c (e, l, j)g c eJ (r s )] , (9) 

at any radius where V(r) has become well represented 
by —Ci/r — C 4 /r 4 . Compared to single-scale QDT for¬ 
mulations [261 . the formula for bound spectrum in the 
two-scale formulation, Eq. ©, is structurally the same 
except that the two-scale depends parametrically on 
Pi/Pi, and the two-scale K c is defined in reference to 
—Ci/r — C 4 /r A , instead of —C\/r solutions. 

From our analytic solutions, we obtain 


tan 9 l +tan(7a//2)(l + M;)/(l-M;) 

1 — tan Op tan(7w/2)(l + M;)/( 1 — Mi) 


Here v is the characteristic exponent for the —l/r a — 
(P 4 /Pi ) 2 /r 4 potential, given in the supplemental mate¬ 
rial [28|. It is a function of the scaled energy e s and 
depends parametrically on P 4 /Pi, as are the tanQp and 
Mi functions in Eq. m- The function tan 6 l is given 
by tan Of = Yf/Xf, with Xf = E^L^-l)-^, and 
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Here b^s are the coefficients 
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for the generalized Neumann expansion 19|, !29j] of 
function. They are given by 6/" = 1, and 
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FIG. 1. (Color online). The fi c representation of the multi¬ 
scale QDT bound spectrum, as the crossing points between 
\i(e 3 ,Pa/Pi) (the curves plotted) and the function fj, c (e,l,j ) 
(the nearly horizontal lines), illustrated here for 2 S' 1 / 2 (1 = 0) 
series of Cs (Black) and Rb (Blue). When the spectrum is 
known, the function evaluated at bound state energies 
gives a discrete representation of the u? function. The exper¬ 
imental data for Cs [jf] and Rb fill . [30] (squares) show that 
and therefore iv c , are to a very good approximation a 
constant, including those for the ground states (6s for Cs and 
5s for Rb). The as a function energy is different for Cs 
and Rb due to their different length scale ratios Pa/P\, which 
in turn is due primarily to their different core polarizabilities. 
a C ore is taken here to be 15.77 a.u. Q for Cs, and 9.076 a.u. 
for Rb [©]. 


/i c (e, l, j), where \i := [tan _1 (xf) — 7 t/4]/7t and is the 
“new” quantum defect defined by := [tan^ATQ — 
7 t/4]/7t. Both Xi and fi c are taken to be within a 
range of [ 0 , 1 ) by taking tan - 1 (x) to be within a range 
of [7 t/4, 57t/ 4). In applying Eq. © to the spectra of 
an alkali atom, one has C\ = 1 and C 4 = a core /2 in 
atomic units, where a C ore is the polarizability of the core 
(i.e. the ionic core excluding the outer electron), and 
Pa/P i = 2^ 3/ V C^core- 

The key difference between the two-scale QDT and the 
standard single-scale QDT i! lies in the fact that the 
“new” quantum defect fj, c , being defined in reference to 
the —C\/r — C^/r^ potential, is determined by the loga¬ 
rithmic derivative of the wave function at much shorter 
distances than the standard quantum defect /Uy oul de¬ 
fined in reference to the —C\/r Coulomb potential, and 
has therefore much weaker energy dependence. Figure [1] 
illustrates both the /z c representation [25| of the two- 
scale QDT spectrum, and the weak energy dependence 
of /z c (e, Z,j), using experimental data for the 2 Si /2 series 


of both Cs Q| and Rb DIE!. 


From a different perspective, the two-scale QDT pro¬ 
vides an analytic description of the energy dependence 
of the /r^ oul in the Rydberg-Ritz formula. Specifically, 
Eq. © for the spectrum can be solved, more precisely 
re-casted, as the solutions of 
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in which to is a positive integer, v 0 := l + 1/2, A := 
(Pi/Pi)y/-e s , w~ := -1/2+ l/(2y/—e s ), and c m (p) := 
11 /lo' Q(. u + j) with Q(v) given by a continued fraction 


QM = 


1 


1 _ ( P±\2 _ ++ 3 / 2 ) 2 + 1/4 _/-)/.. 4 - 1 1 

1 ^Pl> ++l)[++l) 2 -^]++2)[++2)2- I ^] < *'+ + D 

_ _ (13) 

The Mi function in Eq. (HOD is given by Mi := [T(z/ - 
w~)/Y(—v — w~)]Mi , where 


Mi := [(pi/p^eslY 

r(l - I/)T(1 - v + r , p)T(l -v-v 0 ) / Ci(-p) \ 

T(1 + jz)T(1 + v + r , o)T(l + v — vo) \Ci(-\-v)) 

(14) 

in which Ci(v) := fl/lo Q( v + /)• 

Equation © gives the two-scale QDT spectrum as the 
cross points between a universal function Xi( e s, Pi/Pi) 
and a short-range K c (e,l,j) function. It can also be for¬ 
mulated in a n c representation [25|, as x/(e s , Pa/P i) = 


Here e n ijs is a scaled bound state energy (hence the 
subscript s) defined by e n i js := ( E nij - E ion )/s^ = 
(E n ij E[ on ) / A Rm i 

tan( 7 r/iy° ul ) = j- t<m[n(v - 1 / 2 )] , (16) 

1 + ti 

with 


r r(i/( 2 w a ) + i/ + i/ 2 ) ~ 

‘ [ t ( 1 /( 2 k s )-v + 1 / 2 ) \ 

K c cos(7ri//2 — Op) + sin(7TJ//2 — Op) 

I\ c cos(nv/2 + Op) — sin(7ri//2 + Op) 

Equation m, which is formally equivalent to Eq. © 
with a different scaling, shows that the spectrum for any 
potential that behaves asymptotically as —Ci/r— C^/r^ 
can be expressed as a Rydberg-Ritz formula with an 
energy-dependent /x^ oul , consistent with the general con¬ 
clusion by Hartree many years ago |2|. The energy depen¬ 
dence is described analytically by Eq. (fl6l) . It is, to the 
best of our knowledge, the first nonperturbative analytic 
description of this energy dependence. 

Figure [2] compares the Coulomb quantum defects rfp° nl 
to the two-scale QDT predictions of Eq. m with a con¬ 
stant K c , using as examples the 2 S \/2 and 2 Pi /2 series 
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FIG. 2. (Color online). Energy variation of the Coulomb 
quantum defects starting from the ground states. Open sym¬ 
bols: data derived from experimental spectra [&, 3, |7l], j3(| 
and Eq. m- Solid dots: QDT predictions from Eq. m 
with constant K c s, Connecting lines: Eq. Hi) which is well 
defined for all negative energies. 


of Rb and Cs. The Rb data are taken from recent preci¬ 
sion measurement and analysis in Refs. SEUEi The 
a C ore is taken to be 9.076 a.u. 0. The dimensionless 
parameter K c is determined to be —1.314 for the 2 Si / 2 
series, and —1.173 for the 2 P \/2 series, by fitting to an 
intermediate-n portion of the spectra. They correspond 
to ji c = 0.4571 for the 2 S\/ 2 series, and //' = 0.4747 for 
the 2 P \/2 series. The Cs spectrum data are from Ref. Q, 
and we have taken its a cor e = 15.77 a.u. [|]. The pa¬ 
rameters K c are determined to be —0.9135 for the 2 Si/ 2 
series, and —0.8942 for the 2 T’i /2 series. They correspond 
to n c = 0.5144 and = 0.5178 for the 2 £ 1/2 and 2 Pi/ 2 
series, respectively. It is remarkable that even with such 
constant K c s (or /As), the two-scale QDT predicts /iy oul 
with an accuracy better than 0.3% for the ground states, 
and the energy with an accuracy better than 1.3% for the 
ground states, and progressively better for excited states. 

The weak energy dependence of K c (or /A), and the 
degree to which a constant K c (or /A) describes the en¬ 
ergy dependence of Aiy oul show that atomic spectrum 
follows the universal behavior as characterized by the 
—Ci/r — C 4 /r 4 solution, not only for high partial wave 
states [lj, |l5| , but also for the S and P states, and not 
only for highly excited states, but also for the first few 
excited states and the ground state. The weak energy 
dependence of I\ c also implies that the probability for 
finding this outer electron in the region where the po¬ 


tential defers substantially from —Ci/r — Ci/r 4 is small 
pH ], and the wave function, including its normalization, 
is accurately given by the analytic —C\/r — Ci/r 4 wave 
function. This combination of spectrum and wave func¬ 
tion both following a broader universal behavior is what 
will lead to the broader universal behaviors in atomic po¬ 
larizability and the Cq coefficient for different atoms and 
different electronic states. 

Through the weak energy dependence of its short- 
range parameters, the two-scale QDT allows the determi¬ 
nation of the Rydberg spectrum from the measurement 
of the first few excited states, and allows the spectral de¬ 
termination of the core polarizability [§7f. Above the ion¬ 
ization threshold, the same theory provides an analytic 
description of electron-ion scattering i&0 over a wide 
range of energies. Higher accuracy on the spectrum and 
other atomic properties, when desired, can be achieved 
by taking into account the weak energy dependence K c 
(or /A) using a standard Taylor expansion (since they 
are analytic functions of energy, unlike u Coul in the pres¬ 
ence of_the polarization potential) [32j. Multichannel 
i!0 and anisotropic generalizations of the theory 
will extend its description to atomic species other than 
group-I atoms. The two-scale QDT can also be used in 
a fully ab initio fashion together with a i?-matrix theory 
@,0 , leading to more efficient and accurate calculations 
with a smaller R-matrix box. 

In conclusion, we have presented a two-scale QDT for 
a Coulomb plus polarization potential, and have used it 
to establish a broader universality in atomic spectrum, 
covering not only the Rydberg states, but also lower ly¬ 
ing states including the ground state. Beyond single¬ 
atom and two-atom properties (such as Cq), the the¬ 
ory will have implications on molecular electronic struc¬ 
ture, and on atomic and molecular interactions and re¬ 
actions. For molecular structure, for instance, it implies 
that the electronic structures, including both ground and 
excited bands, of Cs,/v and Rb^r clusters differ from each 
other only through different polarizabilities of Cs + and 
Rb + ions and a few quantum defects fi c (l,j) for Cs and 
Rb atoms (with N —> 00 corresponding to condensed 
phases). Mathematically, the same — C\/r — Ci/r 4 so¬ 
lution is applicable not only to electron-ion, but also to 
ion-ion interactions. Finally, this first establishment of 
analytic multiscale QDT gives hope that similar theories 
can be developed for other interactions, such as atom- 
atom [34] and ion-atom [25il425i . that are currently 
treated either at only a single scale 01 S 3 , or at mul¬ 
tiscale but only numerically (see, e.g., Refs, (ilil - lssl ]). If 
successful, such developments will have impact on almost 
every aspect of atomic and molecular structure, interac¬ 
tions, and reactions. 
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